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The twisted reduced model of large QCD with two adjoint Wilson fermions is studied numeri- 
cally using the Hybrid Monte Carlo method. This is the one-site model, whose large limit (large 
volume limit) is expected to be conformal or nearly conformal. The symmetric twist boundary 
condition with flux k is used. ^=0 corresponds to periodic boundary conditions. It is shown that 
the quark mass and dependencies of the model with non-vanishing k differ significantly from 
those of the k=Q model. A preliminary result for the string tension calculated at N=2d,9 is pre- 
sented. The string tension seems to vanish as the physical quark mass decreases to zero in a way 
consistent with the theory being governed by an infrared fixed point with = 0.8 1 .2. 



The 30th International Symposium on Lattice Field Theory 
June 24 - 29, 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the atithor(s) under the terms of the Creative Commons Attribution-NonCommercial-ShitreAlike Licence. 



http://pos.sissa.it/ 



Twisted reduction in large N QCD with two adjoint Wilson fermions 



Masanori Okawa 



1. Introduction 

Recently the SU(N) lattice gauge theories with two adjoint Wilson fermions have received 
much attention, since their massless limits are expected to be conformal or neai^ly conformal[|l|]. 
In fact the first two coefficients of the perturbative j8 function expressed in terms of the 't Hooft 
coupling are given by bQ = {n-ANf)/2AK^ and^i = (17 - 16A^/)/192;r'^, with A^y the number of 
adjoint quark flavors. This would imply the existence of an infrared fixed point if 17/16 < A'^/ < 
11/4 irrespective of the value of gauge group N. In particular, the SU(2) theory has been studied by 
many authors, and now there is a consensus that this theory is conformal and the chiral condensate 
anomalous dimension at the infrared fixed point is approximately 7* ~ 0.3 |[l]]. 

The purpose of the present talk is to study the large lattice gauge theory with two adjoint 
Wilson fermions directly by means of the twisted space-time reduced model||^, ^]. This is a one- 
site model given by four SU(N) matrices only. There is a close connection between this model 
with = L^, with the usual SU(N) theory on a space-time lattice. In the lai^ge L limit this model 
should reproduce the corresponding large A^ infinite volume theory. It should be noted that the use 
of the space-time reduced model has a great advantage in the treatment of the dynamical quark 
effects. It is almost impossible to make a simulation with both a large SU(N) internal group and 
large volume V , since the number of the degrees of freedom is too large. On the other hand, the 
number of degrees of freedom is significantly reduced, by a factor V , in the space-time reduced 
model, and the study of the dynamical quark effects is manageable in this theory. 

It is widely advocated [Q, ^, ^ that the introduction of the dynamical adjoint Wilson fermions 
is enough to restore the Z'*(N) symmetry which is needed for the reduction idea to hold[^], and is 
broken[§] in the original Eguchi-Kawai model (EK model) Our results also confirm that, with 
the inclusion of the dynamical adjoint fermions, the Z^(N) symmetry is indeed unbroken for a wide 
range of values of A'^ and bare quark masses. At the same time, however, we also find that this theory 
has a very strong A'^ dependence, making it practically useless for extracting the large N behavior. 
Our proposal is to apply twisted boundary conditions in the formulation of the space-time reduced 
model. For the pure gauge theory, we have recently demonstrated that the twisted reduction (the 
twisted Eguchi Kawai model) works quite well, allowing a determination of the continuum string 
tension!^, 10 1. We will show that the model with dynamical adjoint Wilson fermions and twisted 



boundary conditions has significantly less A'^ dependence than its untwisted counterpart. This has 
allowed us to present a preliminary calculation of the large N string tension using N=2S9. We show 
that the string tension seems to vanish as we decrease the quark mass towards zero, as predicted 
by the existence of an infrared fixed point. The rate of decrease points towai^ds a large anomalous 
dimension 7* = 0.8 ~ 1.2 at this point. 

2. Formulation 

We consider SU(N) gauge group with N = L?, L being a positive integer. Then the action of 
the model is given by 

4 Nf _ 

S = bN Tv{Z^yU^UyUlul) + Y,'^jDw'i'j. (2.1) 

b is the inverse 't Hooft coupling b = \ /g^N. are four SU(N) matrices and ^'j are the fermion 
fields in the adjoint representation of color SU(N) for Nf = 2. The twist tensor 
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2ni 

c^v = exp 1 k— 



(2.2) 



introduces symmetric twisted boundary conditions to the space-time reduced model specified by 
the flux k. The weak coupling analysis demands that k and L are co-prime for the model to work 
well in that limit [^]. The k=0 case corresponds to periodic boundary conditions. The Wilson Dirac 
matrix is given by 



D 



w 



1 



4 



7,)<^ + (l + 7,)<^" 



(2.3) 



U'ff^ are the link variables acting on the adjoint fermions in the color {N,N) representation as 

We have studied the previous model by means of the Hybrid Monte Carlo method. The sim- 
ulation has been done on Hitachi SR16000 super computer with one node having 32-cores IBM 
powerV and peak speed of 980 GFlops/sec. Thanks to the Hitachi system engineers, our codes 
are highly optimized for SRI 6000 and the sustained seed is 300-600 GFlops/sec depending on the 
values of A'^ and K. 

To study the A'^ dependence of the model precisely, we have performed simulations at A'^= 25, 
49, 81, 121, 169, 225, 289, which con^espond to all odd integer values of L from 5 to 17. It is 
known that part of the finite A'^ effects for our twisted model correspond to the finite volume effects 
in the usual lattice theory on an L"^ lattice[^]: Hence, we based our results on Wilson loops W{R, T) 
up to R=T={L - l)/2. For N=2S9 this corresponds to R=T=S. 



3. fc and dependence 



In fig. 1 we show the K dependence of the average of plaquette E =< Z^yU^UyUjiUl > at 
^=0.35. For periodic boundary conditions ^=0, we show the results for A'^=25,49 and 81. As K" 
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Figure 1: k dependence of E for various at b=Q.35 and ^=0 and 1. 
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b=0.35, K=0.12 
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Figure 2: dependence of E for various k at ^7=0. 35 and k'=0. 12. 



increases E decreases, which is quite unnatural since as K increases one usually expects to have 
larger dynamical quark effects which tends to order the system and, thus, to increase the internal 
energy E. Furthermore, there exists an abrupt change in E at intermediate values of K. However, 
the position of the abrupt change depends significantly on the value of A'^. This change is not 
related to Z^(N) symmetry breaking, which was found to remain unbroken for the set of simulation 
pai^ameters presented in this talk. The large A'^ dependence suggests that it is related to some finite 
A'^ artifact. The results for twisted boundary conditions and k=\ are shown for N=49, 81 and 121. 
In this case E increases as we increase K as expected, and the A^ dependence looks rather small 
within the rough scale of the plot. 

A more precise study of the A'^ dependence is given in fig. 2, where E is plotted as a function 
of \/N for fcO.35, K'=0.12 and various values of k. As we have stated previously, at JC^O. 12 there 
is no Z'*(N) symmetry breaking. However, the results at finite value of K seem to be strongly 
affected by the nearby presence of the Z^(N) symmetry breaking point at k"=0, namely the reduced 
model without fermions. We know that for the pure gauge theory (k"=0) the breaking occurs as 
follows||8],[n|,|l|,[T|,|3l]: 

• ^=0, Z'*(N) symmetry is broken for all N >0. 

• k=l, Z^(N) symmetry is broken for A'^ > 100. 

• ^=3, Z^(N) symmetry is broken for A'^ > 784. 

• ^=5, Z'^(N) symmetry is expected to be broken for A'^ > 2200. 

Now let us look at fig. 2 more carefully. The data with ^=0 have, again, a large A'^ dependence. 
For A'^ <121, this dependence seems linear in 1/A'^. However, such a linear extrapolation to the 
large A'^ value of E is quite misleading, since the behavior of E suddenly changes for A'^>121. The 
data with k=l have less A'^ dependence. However, as we increase beyond A'^=121, E also starts to 
decrease. We think this behavior is closely related to the Z^(N) symmetry breaking at K=0. If this 
interpretation is correct, we should expect a similar pattern for ^=3(5), with a smooth A'^ dependence 
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up to < 1%A{N < 2200). In fact, we observe that the two data points at N=289 for k=3 and 5 ai'e 
quite consistent. In fig. 3 and 4, we display the same type of information but for (b, k) values of 
(0.36, 0.12) and (0.35, 0.14) respectively. Details are different, but the overall trends are the same 
as in Fig. 2. 
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Figure 3: dependence of E at b=0.36 and fc=0.12. Figure 4: dependence of E at ^7=0.35 and K=0.14. 



4. String tension at A'^=289. 

Based on the results of the previous section, we made an attempt at a determination of the 
string tension for N=2S9 and k=5 and various values of K. We used the same procedure as used 
for the pure gauge reduced model in order to extract the string tension from rectangular R ^ T 
Wilson loops W{R,T) ^ |l0|]. Namely, we extract the string tension from Creutz ratios formed out 
of Wilson loops with smeared links. The Ape smeared links, defined as 



jjsrneared ^ p^^-^ 



[/^ + C £ {Zy^UyU^,Ul+Z^yUlU^,Uy 



(4.1) 



are used to enhance the signal to noise ratio. Proj^ stands for the operator for projection onto the 
SU(N) matrices. 

In fig. 5, we show Wilson loops W{R,T) after 5 smearing with smearing parameter c=0.1. We 
observe a clear signal up to loop size 8, which according to our previous comment, corresponds to 
slighly less than half the effective linear lattice size L=17 for our value of N=2S9=n^ . 

In fig. 6, we compare Wilson loops calculated with different k values. In this rough scale, 
Wilson loops with R=l, W{R = l,T), look consistent for all values of ^=0,1,3 and 5. On the other 
hand, the situation is quite different for R=5. The results with k=3 and 5 coincide with each other, 
giving what we think are the correct values. It is amazing that the results with ^=0 are so different 
from those of k=3 and 5 as T increases. The Wilson loops W{R = 5,T) with k=l also have smaller 
values than those with ^=3 and 5. We believe that this is also related to the instabilities occurring 
for A: = 1 at > 100 with K-=0. 

In fig. 7, we show with red symbols our preliminary results for the extracted string tension a 
as a function of K calculated at N=2S9, k=5, Zj=0.35. The value at k:=0 is obtained from the Twisted 
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SU{289), k=5, b=0.35, )t=0.12, after 5 smearing 
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Figure 5: Wilson loop after 5 smearing. 



Figure 6: Wilson loop for various values of k. 



string tension and iowest eigenvalue of Q , SU(2a9), k=5, b=0.35 



string tension and lowest eigenvalue of Q , SU{289), k=5, b=0.36 
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Figure 7: a and lowest eigenvalue of Q ?Xb = 0.35. Figure 8: O and lowest eigenvalue of g at = 0.36. 



Eguchi Kawai model without fermions. As we increase K, the string tension decreases and seems 
to vanish around k:=0. 17. 

Although at present we have not calculated any fermionic spectrum, it is straightforward to cal- 
culate the lowest eigenvalue of the positive square hermitian Wilson-Dirac operator = {ysD-wY, 
which should be related to quark mass square. In fig. 7, we also plot the lowest eigenvalue of 
with green symbols. It seems that the string tension and the lowest eigenvalue of vanish at the 
same point. 

If the massless theory is governed by the infrared fixed point at which the theory is conformal, 
all physical quantities having dimension of mass square should behave as 



1 1 \ W 
a ~ i 

K Kr 



(4.2) 



with the chiral condensate anomalous dimension at the infrared fixed point. We assume that the 
quark mass is proportional to (1/k: — 1/ Kc). Fitting the b = 0.35 data in the range K"=0.12 - 0.16 
with a common vanishing point K" = for the string tension and the lowest eigenvalue of Q^, we 
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get a rather lai^ge estimate for = 0.81(8). This is considerably larger than the value obtained by 
other authors for A'^=2[|]]. It should be noted, however, that our data of the string tension are quite 
preliminary. In particular, it is based on small number of configurations of order 100 with a rather 
low global acceptance ratio in HMC at K"=0.15 and 0.16. 

In fig. 8, we repeat the analysis for the data a.tb = 0.36. Again, requiring that the string tension 
and the lowest eigenvalue of vanish at the same point, our data for the string tension in the range 
JC=0.I2 - 0.15 give y* = 1.17(21). 

5. Conclusions 

We have demonstrated that the twisted reduced model of large A'^ QCD with two adjoint Wilson 
fermions works quite well. The N dependence of the model with twisted boundary conditions 
(k 7^0) is significantly smaller than that of the model with periodic boundary condition (^=0). The 
string tension calculated at N=2i9, k=5, clearly decreases as we increase K and seems to vanish 
around K"=0.17 in a way consistent with the expectations for an infrared fixed point with a rather 
large value for 7* = 0.8 ~ 1.2. 

We have started new simulations for small quark mass with reasonable global acceptance ratio. 
We have also launched the simulation of the single flavor adjoint Wilson fermion using the Rational 
Hybrid Monte Carlo method. We hope to present new results at the coming lattice conference. 
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